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The neutrinoless double beta decay (0^/3/3) induced by light Majo- 
rana neutrino exchange between decaying nucleons, accompanied by the 
squark exchange inside one nucleon, recently discussed by Babu and Mo- 
hapatra, is carefully analyzed both from the particle and nuclear physics 
| sides. New nuclear matrix elements relevant to this mechanism are cal- 

CNJ ■ culated. We extend the analysis to include mixing of light neutrinos 

with heavy and "sterile" neutrinos. It introduces another supersym- 
^ | metric (SUSY) contribution to 0vf3(3. We discuss constraints on the 

0^ ' Ijtp MSSM parameters imposed by the current experimental limit on 

Ov/3/3 decay halfdife of 76 Ge. 

o Neutrinoless double beta (GV/3/3) decay is a sensitive probe of physics be- 

yond the standard model, since it violates lepton number. Recent experi- 
mental progress has pushed the existing half-life limits of decay beyond 
Ti/2(0v(3(3) > 7.4 x 10 24 years and further progress can be expected in the near 
future [|IJ . This experimental result casts stringent constraints on new physics. 
Particularly, for the conventional mechanism of GV/?/3-decay with massive Ma- 
jorana neutrino exchange between decaying nucleons (see fig. 1) it implies an 
upper bound on the neutrino mass below 1 eV 0. (There exist, however, other 
mechanisms which might induce Qv(3(3 decays as well M, |3|].) 

In this paper we study contributions to decay within the R-parity 

violating Minimal Supersymmetric Standard Model (^yVlSSM). The ^yVlSSM 
has been extensively discussed in the literature since it has very interesting 
phenomenological Q and cosmological || implications. It also gives a very 
natural framework for rare lepton number violating processes and particularly 
Ou/3/3 decay 1-0. 



The supersymmetric mechanism of Qvj3(3 decay was first proposed by Mo- 
hapatra and later studied in more details in Refs. 0, ||. In Ref. @ it was 
shown that the gluino exchange contribution to Oz//3/3-decay leads to a very 
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stringent limit on the first generation ftp Yukawa coupling A' m < 3.9 • 10~ 4 . 
Recently, Babu and Mohapatra [10] found another contribution comparable in 
size with the gluino exchange. It allows one to set stringent limits on combi- 
nations of the intergeneration $p -Yukawa couplings such as A' 1H A / lil , where i 
denotes generations. 

However, in ref. [|IJ limits were deduced using a simplified estimation of 



the nuclear structure matrix elements. The approach used in |T(J is based on 
a simple replacement of the virtual particle momenta and energies by the (es- 
timated) Fermi momentum pp and energy Ep of a nucleon inside the nucleus, 
therefore neglecting essentially all nuclear structure effects. Since the stringent 
constraints found in [II| may have important consequences for the _^ P MSSM 
phenomenology, it is desirable to substantiate these results by detailed calcu- 
lations. 

In this letter we present the results of detailed particle and nuclear physics 
calculations relevant to the contribution of the diagrams in fig. 2 to Qv(3(3 
decay. We extend the previous work to include mixing of the standard light 
neutrinos with exotic heavy neutrinos, neutralinos or some other neutral heavy 
particles or light "sterile" neutrino singlets |ll| . Inclusion of mixing leads to 
the new slepton exchange diagram in fig. (2.b). 

The .^yVISSM is an extension of the MSSM which results from the inclusion 
of explicit i?-parity violating (_^ p ) terms W^ p into the superpotential W = 
Wmssm + Wft p , where 



w % = A, ; ,/.,/. ; /-;, • y^L.Qji), • Ay,/;,/;, (i) 

Indices i,j, k stand for generations. L, Q denote lepton and quark doublet 
superfields and E, U, D lepton and up, down quark singlet superfields. The 
first two terms in eq. (|I|) lead to lepton number violation, while the last one 
violates baryon number. For Ov/3/3 decay only the A and A' type couplings are 
of relevance. 

Mixing between scalar superpartners of the left and right-handed 
fermions Jl,r will play a crucial role in our subsequent consideration. It occurs 
due to non-diagonality of the mass matrix which can be written as 

/ m 2 h +m}-0A2D z -m f (A f + fi tan/3) \ 
f " ^ -m f (A f + fi tan/5) m 2 f + m) - 0.08D Z )' [ > 

Here, f = d,s, b, e, fi, t and / are their superpartners. Dz = Mf cos2/3 with 
tan/? = (H®) / (H®) being the ratio of vacuum expectation values of the two 
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Higgs doublets, m; s are soft sfermion masses, Aj are soft SUSY breaking 
parameters describing the strength of trilinear scalar interactions, and \x is the 
super symmetric Higgs(ino) mass parameter. Once sfermion mixing is included, 
the current eigenstates fa, j# become superpositions of the mass eigenstates 
fi with the masses mj and the corresponding mixing angle 9? is defined as 



in 
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m LL + m RR =F \l ( m LL ~ m RR) 
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where m 2 LR , m 2 LL , m 2 RR denote the (1, 2), (1, 1), (2, 2) entries of the mass matrix 
(©)■ 

Now it is straightforward to find the effective 4-fermion v — u — d — e vertex 
induced by the sfermion exchange in the diagrams presented in fig. 2. The 
corresponding effective Lagrangian, after a Fiertz rearrangement, takes the 
form 



eff 
SUSY 



X 



G p 
2rj nj 



I {vUut ~ *n$Ln) ■ Ki ■ (^(1 + 7 5 )e, c ) («(1 + 7b)d) 



+ 
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{l)LL ■ U ni ■ - j 5 )ef) (u(l + 75 )d) + 
\\^rr ■ U m (vi 7 M (1 + 7sh c ) (u 7m(1 - T5)d) + 
\vZ LR ■ U* ni ■ (p, a^(l + 75)S C ) (u <v(l + l5 )d) 



(4) 



The lji p MSSM parameters 77 and neutrino mixing matrix C/y are defined as 
follows 
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YFm- ( 9 ) 



Here f](f)LR denotes the contribution vanishing in the absence of fi — fn. - mix- 
ing while ?7(/)ll and r/^ RR in this limit correspond to the fi and fn exchange 
contribution in fig. 2. We use the notations = d,s,b and e^) = e, fi, r. 
Due to the antisymmetry of the Yukawa coupling \ n jk in nj it follows that 
V^lr = 0- This is an essential difference between the slepton l L — l R and 
the squark — q R contributions. The latter is not imposed to vanish at any 
combination of indexes. 

The matrix element of the SUSY accompanied neutrino exchange mech- 
anism can then be calculated according to the the diagrams in fig. 2 with 
the point-like 4-fermion vertex described by the effective Lagrangian eq. (§) 
with the sfermion exchange parts in the top. The bottom parts of these dia- 
grams is the standard model charged current (SMCC) interaction. Applying 
the standard procedure (for details see ||), one can get the matrix element 
7?-Oi//3^(0 + — > + ) of the Qvf3f3 decay for + — > + transitions. For two outgoing 
electrons in S-wave states it takes the form 

n 0ul3 p(0 + -> 0+) = ^/l[a-e(l + 75)e c + 6-e 7 o75e c ], (10) 
a = (4vt ) LR-V(« )L R)u: i UeM?(m e Ry 1 + (11) 

Trip 



b = ±rffaU«uJ^Mf. (12) 

The term proportional to S ne corresponds to the ordinary neutrino exchange 
contribution with two standard model charged current vertices in fig. 1. The 
normalization factor Cq v is defined as Cq v = (Gp2m e ) /(8y/2irR). We would 
like to stress that the terms proportional to b in eq. (TTD|) can, in principle, be 
discriminated from the terms proportional to a (particularly from the ordinary 
mass mechanism of 0i//9/9-decay) by measuring the angular correlation between 
the two outgoing electrons. 

The following nuclear matrix elements are involved in the calculation of 
Rovpp in eq. ([TUp 



Mf = a x 



M$ + \m%, 



(13) 

Mf = a 2 Mf-M { $ T , M$ = a 3 Mf. (14) 
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They are defined by (summation over nucleons is suppressed) 

(15) 
(16) 
(17) 

0+ > .(18) 

Neutrino potentials can be written in the integral form 

h+(r ab ,m Ui ) 
h R (r ab ,m Ui ) 
h T '{r abl m Ui ) 



Here, ui = yq 2 + rnl.; jk(qr) are spherical Bessel functions and Ro is the 
nuclear radius, introduced to make the matrix elements dimensionless. The 
following notations are used: r ab = (V - ~r*&), r ab = |r o6 |, r ab = r ab /r ab . The 
above formulae have been written in the closure approximation which is well 
motivated for the tiv(3(3 decay 0. A in eq. ( |I9|) is the average intermediate 
state energy. 

Note that the nuclear matrix element Aij), has never been considered in 
the literature before. We have calculated the nuclear matrix elements relevant 
for our subsequent numerical analysis using the pn-QRPA model of [[□[]. We 
take into account both short-range correlations and finite nucleon size effects. 
The latter is described by introducing nucleon form factors in momentum 
space. In the present case these are the SMCC form factors Ffif A) (q 2 ) and the 
isovector scalar and pseudoscalar current form factors Fp S (q 2 ). For all form 
factors we take, as usual, a dipole parameterization: F°-(q 2 )/ Ff(0) = f(q 2 ) = 
(1 + q 2 /m 2 A ) with rriA = 0.85GeV. The q-dependent factor f(q 2 ) is included 
in the definition of the neutrino potentials in eqs. (|19|). The form factor 
normalizations are Fy^(0) = fv(A) ~ 1(1.261). Fg P (0) can be calculated 
within the conventional non-relativistic quark model or the bag model. We 
take their numerical values from ref. [[H -Ff(O) rs 0.48 Fp(0) « 4.41. These 
values correspond to the bag model calculations. 

5 



M 



GT 



M 



GT' 



< 0j\\h + {r ab ,m v .)T+T b + \\0t >, 

< Ql f \\h + (r ab} m Vi )(a a a b )T+T b + \\Q J r >, 

< tf\\h R (r ab ,m Vi )(o a a b )T+T b + \\tt~l >, 

<0f\\h T '(r ab ,m Vi ) (a a r ab )(a b r ab 



T' T h 



2 r°° 
—R I dq- q 

7T Jo 

roo 

/ dq-q 

Jo 



2 

R 

2 R 2 



■Jo(qr a b)f (q 



7T mp 

2 R 2 



TT nip Jo 



dq ■ q 



lu(lu + A) ' 

Jo(qrgb)f 2 (q 2 ) 
oj(u + A) ' 
Jo(qr a b) ~ 3ji(qr ab ) 



u(u + A) 



f{q 2 



(19) 
(20) 
(21) 



The nucleon structure coefficients in (|13|) are defined as a\ = {F P / (2f a), 

« 2 = (fv/fA) 2 , a 3 = (fv/f A )(Fj? ] /f A ). 

Now having the Oz//3/3 matrix element completely specified the inverse half 
life Ti/ 2 (Oi//5/5) can be written as 

TyXOvPP) = \a\ 2 G 01 + \b\ 2 h 9 G 09 + Re(a*b)h 6 G e, (22) 

where he = m e R/8, hg = {m e R) 2 j '16, and phase space factors G^ are given in 
ref. f. 

The formulae presented above describe the contribution of the diagrams in 
fig. 2 for a general neutrino content (see eq. (§)). We now turn to a particular 
case and assume that all neutrino mass eigenstates fall into two groups: light 
neutrinos m Vi < lOMeV and heavy neutrinos m Ui > lOGeV. We denote them 
as Vi and iVj, respectively. There might be also a "sterile" neutrino v s with 
respect to the SM gauge group. Assume further, that the light neutrinos 
have non-negligible mixing with some heavy neutrinos or probably with 
the "sterile" neutrino v s . Than the neutrino composition ([5J) and the unitarity 
relation for the mixing matrix can be written as 

E U V U 3 + Ui 3 U 3 + UisVs, (23) 
3 3 

XX-^ + A,, l±* = Y l U; j U ei + U i8 v a . (24) 

3 3 

Such a structure of the neutrino sector breaks contributions of each nuclear 
matrix elements in eq. ( JTDD into two pieces with a simple dependence on the 
neutrino mass. As an example of this effect consider the first term in eq. (|TTD 

E KiUeiM? = t U* m U ei M\ + t ^ M i- ( 25 ) 

i i i m Ni 

The "sterile" neutrino does not contribute being a SM singlet. The matrix 
elements are M" = M.f{m v = 0) and Mf = lim^^ M N Mf{M N ). They 
do not depend on neutrino masses. Proceeding in a similar way with the other 
matrix elements in eq. ( ]TTj| ) and substituting the result in the half-life formula 
eq. (^) we obtain a polynomial in m v and M^. Assume light neutrinos v to 
be very light while heavy neutrinos iV to be very heavy so that one can neglect 
terms depending on these masses. Than, keeping only leading terms we retain 

TyKOvm = G 01 Ml (m.R)- 1 (4% - %) + V( q )f , (26) 
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where eq. and the property rf^ LR = have been used. We denoted 

V(q) = V(q)LR an d introduced the effective parameters as 

V(l, q )=J2 A nV(l] q) LR (27) 
n 

For the f]^ summation starts from n = 2. The nuclear matrix element hA\ in 
eq. (BB) can be directly obtained from eq. flT3) for A4i as explained after eq. 



(|25|). This matrix element has never been calculated in the literature before. 
Its value calculated in the pn-QRPA for the particular case of 76 Ge analyzed 
below is M[ u \ 76 Ge) = 2.1. [] 

Now we are ready to discuss constraints on the parameters in eq. (|26|) 
imposed by the current experimental lower half-life limit. We use the result 
from the Heidelberg-Moscow 76 Ge experiment Ty^ \ 76 Ge, + — > + ) > 7.4 x 
10 24 years 90% c.l. 

Disregarding a situation with unnatural fine-tuning between the differ- 
ent terms in eq. ( ^6l) one can extract individual limits, numerically tjm < 
2.9 x 10 -8 , f](i) < 7.2 x 1CT 9 . These constraints lead, in principle, to a multi- 
dimensional exclusion curve. A simplified picture can be obtained under some 
reasonable assumptions. Assume all the MSSM mass parameters in eqs. ([5]) 
and (0), (Q) to be approximately equal to the "effective" SUSY breaking scale 
Asusy- Than we get a simplified set of constraints 

^<») 3 ' ^nA^(^) 3 . (28) 

In the last equation we kept only the term corresponding to f exchange. For 
the known values of SM quark and lepton masses in the following we use 
ra^ = rrid = 7.5MeV, = m s = 150MeV, rrid 3 = nib = 4.5GeV and 
m T = 1.8GeV. These masses are present in eqs. (|2|), ©. Then, e' 123 = 
{6.4 x KT 5 , 3.2 x KT 6 , 1.1 x 1(T 7 } and e = 6.5 x 1(T 8 . 

To obtain information on the ^Yukawa couplings themselves one may use 
some reasonable values for Asusy- If one takes, following ref. flUf] , the value 
Asusy = 100 GeV lying in the region of the current experimental lower bound 
for the SUSY particles, one gets A' 113 A' 131 < 1.1 x 10~ 7 , A' 112 A' 121 < 3.2 x lO" 6 , 
A'ln < 6.4 x 10~ 5 and A n A'A < 6.5 x 10 -8 . The former two limits correspond 
to those in ref. JTTJ, but are somewhat weaker than quoted by Babu and 



Mohapatra. This difference can be partly traced back to our numerical value 



1 Recall that this numerical value corresponds to our dimensionless convention. 
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of the matrix element, which turns out to be smaller than anticipated. The 
limit on X'1 n can be compared with the corresponding limit obtained from 
the gluino exchange diagram in ref. [^]. The latter being \'j n < 1.52 • 10- 7 is 
much more stringent then the one derived here. More conservative estimations 
can be derived from eq. ( f28l) implying Asusy ~ ITeV motivated by the SUSY 
naturalness arguments. Then, A' 113 A' 131 < 1.1 x 10~ 4 , A' 112 A' 121 < 3.2 x 1CT 3 , 
A 'in < 6.4 x 10- 2 and A n A'A < 6.5 x 10" 4 . 

Neglecting mixing between the light SM non-singlet neutrinos and the non- 
standard sector discussed above, one arrives at the case considered by Babu 
and Mohapatra [IU|. In this case A n = and in turn fj(l,q) = 0. Then the 



squark exchange diagram in fig. 2(a) is the only contribution to eq. (p6|). 
Introducing mixing one can obtain new information about A type interactions. 
In this case, however, this information is accessible only in a form of some 
effective value of the l/l p Yukawa couplings A A' A. The upper limit for the 
couplings itself implies an extra uncertainty due to the unknown mixing factor 
A. 

The new contributions analysed in this letter together with the previously 
discussed in the literature ||-f| complete the tree- level mechanisms of Ovflfl 
decay within the R-parity violating Minimal Supersymmetric Standard Model 
(-^pMSSM). A detailed presentation of the analysis will be given elsewhere. 

In conclusion we would like to stress that the current experimental limit 
on Oz//3/3 decay half life allows one to establish rather stringent limits on 
l/lp violating SUSY interactions. This might be an additional valuable mo- 
tivation for present and forthcoming Ozy/3/5 decay experiments. 
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Figure Captions 

Fig.l Feynman graphs for the conventional mechanism of 0vf3(3 decay by ex- 
change of a massive Majorana neutrino. 

Fig. 2 Feynman graphs for the supersymmetry accompanied Majorana neutrino 
exchange mechanism of the Ovflfl decay corresponding to (a) squarks and 
(b) sleptons contributions. 
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Figure 2. a 




